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*A  THEORY  FOR  THE  LAMINAE  SUBLAYER 
OP  A  TURBULENT  FLOW 

ABSTRACT 

The  so-called  laminar  sublayer  is  shown  to  be  the  region  where  the 
turbiilent  velocity  fluctuations  are  directly  dissipated  by  viscosity.  A 
simplified  linearized  fom  of  the  equations  of  motion  for  the  turbiilent 
fluctuations  is  used  to  describe  the  turbulent  field  between  the  wall 
and  the  fully  turbulent  part  of  the  flow.  The  mean  flow  in  the  sublayer 
and  the  turbulence  field  outside  the  sublayer  are  assijmed  to  be  known 
from  the  experiments.  The  thickness  of  the  sublayer  arises  naturally  in 
the  theory  and  is  directly  analogous  to  the  inner  viscous  region  for  the 
fluctuations  in  a  laminar  flow.  It  is  shown  that  the  large  scale  fluctu¬ 
ations  containing  most  of  the  turbulent  energy  are  convected  downstream 
with  a  velocity  characteristic  of  the  middle  of  the  boundary  layer..  Thus 
Taylor’s  hypothesis  does  not  apply  to  these  large  scale  fluctuations  near 
the  wall.  The  convective  velocity  found  in  the  measurements  of  pressure 
fluctuations  at  the  boundaries  of  turbiilent  flows  is  in  accord  with  the 
theory.  Csilculatlons  are  given  for  the  energy  spectra  and  u'  fluctuation 
level  in  the  sublayer  and  other  aspects  of  the  fluctuation  field  are 
discussed.  It  is  shown  that  the  production  of  turbulent  energy  is  a 

♦A  preliminary  account  of  this  work  was  presented  to  the  Annvial  Meeting 
of  the  Fluid  Bynamlcs  Division,  American  Physical.  Society,  in  November 
1959  at  Ann  Arbor,  Michigan. 
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maximm  where  the  laminar  shearing  stress  is  equal  to  the  turbulent 
sheajrlng  stress.  The  linear  pressxire  fluctuation  field  at  the  edge 
of  the  sublayer  is  calculated  and  found  to  be  much  larger  than  the 
non-linear  field.  Examining  the  effect  of  strong  free  stream  turbulence 
on  lami  nar  boimdary  layer  transition,  it  appears  that  the  physical  model 
underlying  Taylor's  parameter  is  incorrect. 
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1.  INTRODUCTION 


The  laminajr  sublayer  has  been  a  subject  of  controversy  and  investi¬ 
gation  for  more  than  20  years.  The  reason  for  this  interest  is  that  the 
nature  of  the  flow  close  to  the  wall  has  an  Important  influence  on  the 
heat,  mass,  and  momentum  transfer  from  the  boundary.  Furthermore, 
experiments  have  shown  that  the  flow  of  energy  from  the  mean  flow  to  the 
turb\ilent  motion  is  a  maximum  inside  the  sublayer.  Biis  fact  sviggests 
that  an  understanding  of  the  structure  of  tiirbulence  in  a  shear  flow  may 
depend  on  an  understanding  of  the  flow  near  the  wall. 

The  origlnail  idea,  Taylor  (I916),  was  that  in  a  tm*bulent  flow  there 

ought  to  be  a  thin  fluid  layer  next  to  the  surface  free  of  turbulent 

motion,  a  true  laminar  layer.  Studies  of  the  stability  of  Couette  flow 

(the  flow  between  a  fixed  wall  and  a  moving  wall)  had  shown  that  there 

was  a  critical  Reynolds  number  Re  «  500  below  which  adl  eddies  wo\ild  die 

out.  The  crlticad.  Reynolds  number  Uh/v  was  formed  using  the  velocity  U 

of  the  moving  wall  and  the  sepeiration  distance  h  where  v  is  the  klnenatie 

viscosity.  It  was  postulated  that  the  flow  next  to  the  wall  was  equivalent 

to  a  Couette  flow.  The  laminar  sublayer  thickness  6^  could  then  be 

estimated  by  substituting  6  for  the  separation  distance  in  the  stability 

s 

analysis . 

In  1952,  Page  and  Townend  studied  the  fluctuation  field  close  to 
the  surface  using  an  ultramicroscope  for  following  minute  particles  in 
water.  Uiey  found  no  evidence  of  an  eddy-free  region  neea*  the  wall.  An 
Interesting  discussion  of  this  work  is  given  by  Taylor  (1952)  who  examined 
the  possible  connection  between  special  types  of  disturbances  where  the 
velocity  distributions  were  known  and  the  turbvilent  fluctixation  field 
found  by  Page  amd  Townend.  These  experimental  results  were  confirmed  by 
later  hot-wire  measurements  in  air  by  Laufer  (1950,  1953)  and  KLebanoff 
(1954).  Instead  of  being  eddy  free,  the  txirbulence  level,  as  given  by 
the  ratio  of  u’  the  root-mean-sqviare  value  of  the  velocity  fluctuation  in 
the  flow  direction  to  the  local  mean  velocity  U-,  reached  a  wiRvInnim  value 


5 


of  approximately  .4  close  to  the  wall.  Also  the  turbulent  shear  stress, 
as  deduced  fran  the  mean  flow  measurements,  did  not  vanish  in  a  thin 
region  next  to  the  wall  but  instead  varied  continuously  from  zero  at  the 
wall  to  the  level  of  the  wall  shear.  Biiis  It  has  been  clear  for  sane 
time  that  a  theory  of  the  "laminar"  sublayer  m\ist  account  for  the  fact 
that  the  flow  is  turbulent  all  the  way  to  the  wall. 

There  Is  now  a  relative  wealth  of  ec^terlmental  information  on  the 
fluctuation  field  close  to  the  wall  of  a  turbulent  flow.  What  is  needed 
is  a  theoretical  structure  that  will  provide  a  rational  foimdation  for 
the  \inderstanding  and  interpretation  of  the  experimental  observations. 
Several  recent  attempts  have  been  made  to  develop  phenomenological  models 
for  the  flow  in  the  sublayer.  On  the  basis  of  some  observations  using 
dye  in  water,  Einstein  and  Li  (1956)  were  led  to  postxolate  the  periodic 
growth  and  decay  of  a  true  laminar  region  near  the  wall.  An  equivalent 
model  has  been  proposed  by  Haniatty  (1956).  However,  in  all  these  cases, 
agreement  with  the  measurements  of  the  mean  cuad  fluctuation  field  is 
sensitive  to  the  choice  of  critical  parameters  as  well  as  to  arbitrary 
and  scaaetlmes  inconsistent  assumptions  concerning  the  physical  processes. 
The  pxupose  of  this  paper  is  to  make  a  start  toward  the  development  of 
a  theory  for  the  sublayer  which  follows  from  the  Navler-Stokes  equations 
without  the  need  for  phenanenologlcal  assumptions  or  tinkering  with 
adjustable  peirameters. 

We  have  already  suggested  that  the  sublayer  is  only  a  special  part 
of  the  general  turbulent  fluctuation  field  in  a  bounded  shear  flow.  The 
aim  of  the  theory  will  be  to  describe  as  far  as  possible  the  direct 
influence  of  the  wall  on  the  fluctuation  field.  The  txjrbulent  field  out¬ 
side  of  this  region  of  direct  influence  is  assumed  to  be  known  on  the 
basis  of  the  experimental,  measurements.  This  knowledge  is  essential  in 
the  development  of  the  sublayer  analysis.  Fortunately  we  need  not  be 

concerned  with  shear  flow  turbulence  in  all  its  complexity,  since  only 

\ 

certain,  fairly  simple,  aspects  of  such  flows  are  significant  for  this 
problem.  However,  some  of  the  simple  features  have  been  obscured  by 
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the  accepted  methods  for  presenting  the  experimental  results.  It  is 
customary  to  detenaine  the  scales  or  vave  numbers  for  the  turbulent 

measurements  by  using  Taylor's  hypothesis  (Taylor  1958)  to  justify  the 
necessary  space-time  transformation.  Our  first  step  will  be  to  show  that 
this  hypothesis,  which  was  introduced  to  represent  the  turbulence  behind 
a  grid  In  a  wind  t\mnel,  Is  not  valid  in  a  shear  flow  especially  near  a 
wall.  This  euoalysls  will  provide  a  basis  for  re-interpreting  the  experi¬ 
mental  measurements  6ind  will  make  possible  some  important  simplifications 
in  the  subsequent  development  of  the  theory. 
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2.  TAYLOR'S  HYPOTHESIS  Hff  A  SHEAR  PLOW 


Taylor's  hypothesis,  which  has  been  amply  verified  for  a  uniform 
low  turbvilence  flow  involves  two  assumptions: 

a)  The  turbxolence  pattern  is  convected  past  the  measuring  point 
with  the  local  mean  speed. 

b)  The  tiirbulent  fluctviatlng  velocities  are  smal  1  enough  compeured 
to  the  mean  motion  to  insure  little  change  in  the  shape  of  an  eddy  as 
it  is  carried  past  a  fixed  point. 

The  xxse  of  this  hypothesis  in  a  shear  flow  has  previously  been 
questioned  by  Lin  (1955).  EssentieLLly,  Lin  investigated  the  conditions 
for  negligible  eddy  distortion,  and  showed  that  there  is  "no  general. 
Justification  of  extending  Taylor's  hypothesis  to  the  case  of  shear 
flow".  He  foxmd  that  unless  an  eddy  component  had  a  scale  much  less 
than  the  boundary  layer  thickness,  it  would  suffer  significant  distortion 
due  to  the  mean  flow  shearing  motion  while  being  carried  past  the  measuring 
point  by  the  mean  flow.  However,  Lin's  analysis  did  not  lead  to  an 
alternate  procedure  for  determining  the  turbulence  scales. 

We  will  show  that,  in  general,  assumption  (a)  cannot  be  valid  in  a 
turbulent  shear  flow.  The  departure  from  this  assumption  in  a  boundary 
layer  is  especially  significant  near  the  wall. 

Consider  a  txirbulent  bovuadary  layer  on  a  flat  plate.  At  any  instant 
t,  the  turbulent  fluctuation  field  in  a  boundary  layer  can  be  represented 
by  a  distribution  of  disturbance  vorticity  caaponents  I,  tj,  and  ^ 
throughout  the  boundary  layer.  At  the  wall  the  vertical  perturbation 
velocity  v  must  vanish.  This  boundary  condition  can  be  satisfied  by 
adding  an  image  vorticity  distribution  on  the  opposite  side  of  "Uie  wall. 
Now,  associated  with  the  vorticity  at  a  point  P'  in  the  boundary  layer, 
there  is  an  induced  velocity  at  the  point  P.  The  total  velocity  pertur¬ 
bation  at  point  P  at  any  instant  can  then  be  found  by  integrating  over 
the  boundary  layer  and  image  system  perturbation  vorticity  fields.  The 
extent  of  the  region  over  which  the  integration  must  be  carried  out 
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depends  on  the  scale  of  turbulence.  If  only  small  scetle  motions  are 
present  then  the  region  of  integration  cam  be  confined  to  the  vicinity 
of  P.  There  should  not  be  ajiy  significant  correlation  between  velocities 
at  point  P  8ind  random  small  scale  vorticity  at  distances  from  P  many 
times  the  scale  of  the  dist^2rbances .  If  large  scale  motions  are  present 
then  the  integration  must  extend  at  least  over  the  disteuaces  where  these 
large  scale  motions  are  significantly  correlated. 

A  typical  one  dimensional  energy  spectrimi  for  the  velocity  perturbation 
u  in  the  flow  direction  at  y/5  =  .58  near  the  center  of  a  boundary  layer 
is  shown  in  Fig.  1,  where  y  is  the  distance  from  the  wall  and  5  is  the 
boundary  layer  thickness.  If  we  can  assume  for  the  moment  that  this  dis¬ 
turbance  field  is  being  carried  along  by  the  mean  flow  in  accordance  with 
Taylor's  hypothesis,  then  the  frequency  can  be  converted  into  a  measure 
of  the  space  scale  L  by  the  relation  L  =  where  is  the  velocity  of 

the  local  mean  flow  (we  will  show  that  this  is  justified  in  the  central 
region  of  the  boxmdary  layer). 

Figure  1  also  shows  the  contribution  to  the  total  energy  u  as  a 
function  of  the  scale  of  the  motion.  It  is  evident  that  fully  half  of 
the  energy  is  contributed  by  turbulence  whose  scale  is  more  than  twice 
the  boundary  layer  thickness,  as  first  noted  by  Townsend  (1951).  Thus 
the  velocity  perturbation  at  a  point  P  does  not  depend  significantly  on 
the  vorticity  in  the  immediate  vicinity  of  P,  but  rather  on  the  vorticity 
over  an  extensive  region  of  the  boundary  layer. 

Since  the  vorticity  travels  with  the  fluid  particles,  the  apparent 
velocity  with  which  a  disturbance  sweeps  past  the  measuring  point  P  may 
therefore  be  substantially  different  from  the  local  mean  velocity  at  P. 

For  the  large  scale  motions,  this  disturbance  velocity  will  correspond 
to  the  mean  velocity  near  the  middle  of  the  boundary  layer.  Therefore 
for  points  P  close  to  the  wall,  where  the  mean  velocity  is  low,  we  would 
expect  the  disturbance  velocity  to  be  greater  than  the  local  mean  velocity; 
for  points  P  near  the  outer  edge  of  the  boundary  layer,  the  disturbance 
velocity  should  be  less  than  the  local  mean  velocity. 
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The  space-time  correlation  data  of  Favre,  Gaviglio,  and  Dumas 

(1957,  1958)  provide  a  basis  for  determining  whether  these  deductions 

are  correct.  If  F(f )  represents  the  percent  of  turbulent  energy 

associated  with  the  frequency  f,  and  Taylor's  hypothesis  is  satisfied 

then  the  auto- correlation  coefficient  R  of  the  u'  fluctuation  at  the 

points  P  and  (P+x),  can  be  written  as 

00 


F(f)  cos 


f 


df 


where  R  = 

u  (P) 

T3ie  longitudinal  correlation  coefficient  calculated  in  this  way  should 
then  agree  with  the  longitudinal  correlation  coefficient  directly  measured 
with  two  hot  wires. 

This  autocorrelation  coefficient  has  been  calculated  using  the  energy 
spectrum  measured  at  one  of  Favre's  test  points  and  is  shown  in  Figure  2. 
15ie  integration  is  only  carried  out  to  f^^^  =  1000  c.p.s.,  but  this 
includes  98^6  of  the  energy.  In  Figure  2,  two  additional  correlation 
curves  are  show  for  which  the  integration  was  terminated  at  smaller 
values  of  the  frequency  therefore  eliminating  the  contribution  of  the 
small  eddies.  A  frequency  of  i|00  c.p.s.  corresponds  to  a  longitudinal 
scale  of  about  2. If  centimeters  which  comi«.res  with  a  boundary  thickness 
of  centimeters.  It  is  evident  that  the  correlation  coefficient  at 
large  distances  is  not  significantly  affected  by  the  small  eddies.  We 
shall  confine  our  attention  to  the  portion  of  the  correlation  curve 
determined  by  the  large  scale  eddies  so  that  will  be  lees  than  .k. 

If  the  velocity  of  the  disturbance  U  ,  associated  with  these  large 

w 

scale  eddies  differs  from  the  local  mean  velocity  U.,  U  rather  than  U, 

W  £ 

must  be  used  in  the  formula  for  computing  the  correlation  coefficient 
from  the  energy  spectrum  at  a  fixed  point.  Thus,  at  a  fixed  value  of 
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the  computed  correlation  ciorve  should  he  shifted  horizontally  where 
the  new  horizontal  coordinate  x'  =  x.  f  Otily  portion  of  the 
curve  dominated  by  the  large  eddies  should  be  shifted  in  this  way, 
since  as  the  eddy  size  is  reduced,  the  disturbance  velocity  approaches 
the  local  mean  velocity. 

One  of  the  figxares  from  Favre's  paper  (1958)  is  reproduced  in  Figure 
3.  Auto-correlation  curves  using  Taylor's  hypothesis  and  longitudinal 
space  correlations  measured  with  two  hot  wires  are  shown  for  four  positions 
across  a  bovindary  layer.  While  there  is  the  usual  experimental  scatter, 
there  is  a  systematic  difference  between  the  two  sets  of  curves  below 
R^  w  .4  depending  on  the  location  of  the  meas\iring  point  in  the  boundary 
layer.  Close  to  the  wall  the  measured  longitudinal  correlation  curves 
are  to  the  ri^t  of  the  calculated  auto-correlation  curves.  Near  the  outer 
edge  of  the  boundary  layer,  the  measured  longitudinal  correlation  curve 
is  displaced  in  the  opposite  direction.  At  y/5  =  .24  the  difference 
between  the  two  curves  is  lost  in  the  scatter  of  the  data. 

As  far  as  they  go,  these  measurements  are  consistent  with  the  picture 
in  which  the  laxge  scale  disturbances  which  contain  most  of  the  energy 
move  down  stream  at  a  mean  velocity  characteristic  of  the  central  region 
of  the  boundary  layer  fluid.  At  y/&  =  .24  in  Favre's  bomdary  layer, 
the  velocity  U«  is  approximately  equal  to  .78  of  the  free  stream  velocity 
U^.  The  horizontal  shift  in  the  correlation  curve  that  would  be  expected 
at  y/5  =  .05  if  the  disturbance  velocity  were  equal  to  .78  is  also 
shown  in  Figure  2.  This  shift  is  approximately  the  same  magnitude  as 
the  shift  found  in  Figure  3*  A  similar  divergence  of  the  calculated 
auto- correlation  and  measured  longitudinal  correlation  coefficients  at 
large  scales  has  been  observed  by  KLebanoff  and  led  him  to  remark  that 
"this  divergence  gives  rise  to  the  interesting  speculation  that  the 
large  scale  motions  have  their  own  characteristic  velocity  different 
from  the  mean  speed" . 
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We  conclude,  on  the  basis  of  this  euialysis,  that  in  a  boundary- 
layer,  or  in  fact,  in  any  sheeir  flow,  the  disturbance  velocity  at  a 
measuring  point  P  is  in  general  different  from  the  loceU.  mean  velocity. 
Iherefore,  the  customary  conversion  of  experimental,  spectral  measuronents 
Into  wave  nimbers  Is  in-yalid  in  a  boundary  layer  anywhere  near  the  wall 
except  for  the  small  scale  structure  of  the  turbulence. 

In  the  following  sections  of  this  paper,  many  of  the  numerical  results 
will  be  based  on  the  experimental  measurements  of  KLebanoff.  The  limited 
auto- correlation  emd  longitudinal  correlation  measurements  he  made  do  not 
extend  to  large  enough  scales  to  es-tabllsh  a  -value  for  the  disturbance 
velocity  of  the  large  eddies  with  any  precision.  Accordingly,  we  shall 
use  the  generaJ.  information  obtained  from  Pavre's  da-ta  and  somewhat 
arbitrarily  set  for  the  large  eddies.  This  is  the  -value  of  the 

mean  velocity  at  y/6  =  .27.  We  will  also  need  to  es-bablish  an  approximate 
upper  limit  to  the  frequency  range  to  which  this  dlsturbemce  velocity 
applies.  KLebanoff 's  auto  and  logitudinal  correlation  measurements  Indicate 
that  for  R„  <  *5,  the  expected  shift  of  the  longitudinal  correlation  curve 

X 

with  respect  to  the  auto-correlation  curve  will  have  occurred.  Calculations 
for  as  a  fxmction  of  frequency  at  different  separations  are  shown  in 
Figure  U  for  y/6  =  .05.  If  we  consider  the  curve  for  =1.8,  it  is 

evident  that  -the  frequencies  >  500  make  a  minor  contribution  to  the 
correlation  coefficient.  That  is,  the  frequencies  be-tween  500  and  1000 
only  produce  a  moderate  oscillation  about  the  final  -value  of  R^  =  .45.  A 
frequency  of  500  c.p.s.  corresponds  to  an  eddy  scale  of  L  =  =  4  cm  or 

approximately  (l/2)6.  (Frequencies  between  0  <  f  <  500  c.p.s.  account 
for  about  80^  of  the  fluctuation  energy).  We  will  therefore  limit  the 
specification  =  .8u^  to  frequencies  from  0-^500  c.p.s.  For  high 
frequencies,  or  small  sceLLe  motions,  the  disturbances  move  downstream 
with  -the  local  mean  velocity.  A  very  crude  guess  for  the  dependence  of 
disturbance  velocity  on  frequency  at  y/8  =  .05  will  be  given  in  the  section 
on  the  microscales .  But  KLebeinoff '  s  experimental  data  do  not  pro-vide  amy 
clear  basis  for  guessing  at  this  dependence  in  the  sublayer.  As  will  be 
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mentioned  in  the  section  on  the  pressure  field,  it  may  he  possible  to 
establish  the  variation  of  disturbance  velocity  with  frequency  for  the 
hi^*er  frequencies  fran  the  measurements  of  wall  pressure  fluctuations. 
But  at  the  present  this  information  is  not  available  and  this  lack  of 
knowledge  will  necessarily  limit  certain  possible  applications  of  the 
present  theory. 
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5.  THE  FLUCTUATION  FIELD  OUTSIDE  OF  THE  SUBLAYER 


When  the  experimental  measurements  themselves  are  examined,  they 
reveal  a  rather  striking  similarity  of  the  energy  spectra  at  different 
points  in  the  fully  turbulent  part  of  the  flow.  The  normalized  energy 
spectra  across  a  boundary  layer  as  measured  by  KLebanoff  axe  shown  in 
Figure  5.  Over  the  inner  half  of  the  boundary  layer,  in  the  region  free 
of  intermlttency,  the  spectra  for  the  energy  containing  eddies  appear  to 
agree  within  the  experimental  error.  Differences  in  the  high  frequency 
end  of  the  spectra  would  be  revealed  by  using  a  log  scale  rather  than  a 
linear  scale,  but  these  portions  of  the  spectra  provide  a  negligible 
contribution  to  the  total  fluctuation  energy.  The  same  results  axe 
found  for  pipe  or  chaxinel  flow  spectrail  data.  This  similarity  of  the 
energy  containing  portion  of  the  frequency  spectra  is  just  what  would 
be  expected  from  the  previous  analysis  since  the  disturbeuice  velocity 
for  the  large  scale  eddies  should  not  vary  significantly  across  the  sheax 
flow.  At  each  point  across  the  boundary  layer,  the  hot-wire  probe 
responds  to  disturbances  associated  with  the  same  large  scale  eddy  pattern 
and  so  the  spectral  distribution  should  be  similar. 

On  the  other  hand,  the  fluctuation  energy  varies  frcan  point  to 
point.  Figure  6  shows  the  variation  of  the  root-mean- square  fluctuation 
velocities  u',  v',  and  w'  across  the  boundary  layer.  It  can  be  seen  that 
the  u'  fluctuation  level  increases  by  a  factor  of  2  between  y/B  =  .6  and 
y/6  =  .05.  Why  does  the  fluctuation  level  vary  if  we  are  measuring 
perturbations  due  to  the  same  large  scale  eddy  system? 

A  simple  explanation  for  these  observations  can  be  suggested  by 
considering  the  effect  of  the  wall  on  the  perturbation  field.  If  no  wall 
were  present,  the  induced  velocity  field  associated  with  vorticity  at  P' 
would  be  symmetrical  about  P' .  When  the  wall  is  present,  the  Induced 
velocity  field  of  the  image  vorticity  for  P'  must  be  added  to  the  field 
directly  due  to  P ' .  This  will  cause  an  increase  of  the  induced  velocities 
between  P'  and  the  wall  and  a  decrease  of  the  induced  velocities  beyond  P'. 
Ihis  image  effect  may  be  the  reason  why  the  fluctuation  level  associated 
with  the  same  large  scale  vorticity  Increases  towards  the  weQ.!. 
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So  far  we  have  managed  to  avoid  specifying  the  extent  of  the 
sublayer.  This  has  been  a  somewhat  ambiguous  question  ajid  has  been 
subject  to  different  Interpretation  by  different  authors.  In  Figure 
7  the  variation  of  the  u'  fluctuation  and  the  turbulent  shear  stress 
puv  near  the  wall  are  shown  for  a  boundary  layer  and  a  pipe  flow.  The 
variation  of  puv  has  been  csJ-culated  using  the  measiired  mean  velocity 
profile  and  the  fact  that  the  total  shear  stress  is  essentially  constemt 
near  the  wall.  As  is  custcmary,  the  data  are  presented  in  terms  of  the 
friction  velocity  where  Is  the  friction  at  the  wall  and 

p  is  the  density.  In  both  cases,  the  peak  of  the  u'  fluctuation  is  found 
at  U^.y/v  =  U^.r/v  =  15-20.  Some  authors  specify  U^y/v  «  12  as  the  edge 
of  the  1  flJn^ riftT*  sublayer  and  consider  the  region  12  <  U^.y/v  <  60  to  be 
some  sort  of  transition  region  between  the  laminar  sublayer  and  the  fully 
developed  turbulent  part  of  the  flow.  The  fact  that  the  mean  velocity 
profile  is  nearly  linear  up  to  U^.y/v  12  seems  to  support  this 
definition.  On  the  other  hand,  at  U^.y/v  «  12,  the  turbulent  shear  stress 
is  still  only  about  l/2  of  the  level  of  shear  at  the  wall,  and  only 
asymptotically  approaches  the  wall  shear  value  somewhere  around  U^.y/v  «  100  or 
U^r/v  60  according  to  Figure  7* 

Viewing  the  region  near  the  wall  as  a  whole  we  can  describe  what  is 
observed  in  the  following  general  terns.  Outside  of  the  wall  region,  the 
turbulent  shear  stress  and  the  u'  fluctuation  vary  slowly  compared  to  the 
variations  that  are  found  in  the  wall  region.  Entering  the  wauLl  region 
the  u’  fluctuation  first  rises.  In  Figure  7»  the  fluctuation  level  is 
normalized  in  terms  of  the  fluctuation  level  at  the  edge  of  the  wall  region. 

(it  is  more  difficult  to  identify  the  edge  of  the  wall  region  for  the 
boundary  layer  case.  This  Is  because  the  variation  of  u'  outside  the  wall 
region  is  much  larger  in  the  boundary  layer  case  than  in  the  pipe  flow). 
Apparently,  the  fluctuation  level  increases  mors  in  KLebanoff 's  experiments 
than  in  Laxifer's  case.  However,  the  fluctuation  level  in  the  sublayer 
should  probably  be  compared  with  an  extrapolation  to  the  wall  of  the  u' 
variation  outside  the  sublayers  This  would  make  a  significant  difference 
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in  the  case  of  the  boundary  layer  and  vould  suggest  that  the  rise  In 
fluctuation  level  in  the  boundary  layer  is  not  as  great  as  it  first 
appesLTS  to  be.  Approximately  at  the  point  where  u'  starts  to  rise,  the 
shear  stress  starts  to  decrease  slowly.  The  rapid  decrease  in  u'  is 
confined  to  the  inner  25^  of  the  wall  region.  In  the  theory  given  in 
this  paper,  the  sublayer  is  the  entire  region  between  the  wall  and  the 
fully  developed  turbulent  part  of  the  flow.  There  is  no  theoretical 
distinction  between  an  inner  "laminar*'  sublayer  and  a  transition  region 
although  the  rapid  changes  do  primarily  occur  in  the  inner  portion  of 
the  sublayer. 

The  physical  picture  of  the  large  scale  eddies  containing  most  of  the 
turbulent  energy  moving  downstream  at  a  velocity  of  the  order  of  .8  of 
the  free  stream  velocity  is  reminiscent  of  the  physical  picture  of 
oscillations  in  a  laminar  boundary  layer.  In  that  case  typical  waves 
move  downstream  with  a  velocity  U  (2/5 )U,  and  have  wave  lengths  of 
the  order  of  (2-5)  &.  In  the  equations  of  motion  for  the  perturbations 
in  a  laminar  flow,  the  term  representing  the  action  of  viscosity  is 
negligible  except  in  two  limited  regions  of  the  boundary  layer,  the 
critical  layer  and  the  inner  viscous  layer  close  to  the  wall.  We  shall 
show  that  this  inner  viscous  layer  for  fluctuations  in  a  laminar  bovindary 
layer  corresponds  directly  to  the  sublayer  for  a  turbulent  flow.  Thus 
the  sublayer  is  the  region  where  the  turbulent  fluctuations  in  the  shear 
flow  are  damped  by  viscosity. 

The  existence  of  a  "dissipation  layer"  near  the  wall  was  first 
suggested  by  Townsend  on  the  basis  of  a  study  of  the  turbulent  energy 
balance  in  the  bo\jndary  layer.  He  concluded  that  the  bulk  of  the  turbu¬ 
lent  energy  dissipation  takes  place  by  direct  viscous  action  on  the  large 
eddies  in  a  layer  which  he  thought  was  "most  probably  in  contact  with  the 
laminar  sublayer" .  In  the  present  theory  we  find  that  the  sublayer  itself 
is  a  dissipative  region  and  that  its  structure  is  primarily  determined  by 
the  large  scale  fluctuations  in  the  turbulence. 
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4.  EQUATIONS  FOR  THE  FLUCTUATION  FIELD 


The  equations  of  motion  for  the  fluctuations  in  a  turbulent  field 
(Lin  1959^  P*246)  are  obtained  by  subtracting  the  well  known  Reynolds 
equation  for  the  mean  flow  frm  the  complete  Navier-Stokes  equations. 

We  restrict  our  attention  to  a  steady  flow  in  which  the  mean  velocity 
only  has  a  component  parallel  to  the  flow  direction,  so  that  U  =  U(y) 
and  V  =  W  =  0,  where  U,  V,  and  W  are  the  three  ccmponents  of  the  mean 
motion.  We  will  also  assume  that  the  statistical  properties  of  the 
turbulent  field  such  as  u  6«id  uv  only  vary  with  y.  These  assumptions 
are  reasonably  well  satisfied  by  a  two  dimensional  boundary  layer  flow 
or  a  pipe  flow.  Then  if  u,  v,  and  w  are  the  disturbance  velocities  the 
ccsnponents  of  total  velocity  are  u  =  U+u,  v  =  v,  w  =  w,  and  the  pressure 
isp=P  +  p.  We  can  then  write  down  the  three  equations  of  motion  for 
the  fluctuating  field 


In  each  equation  the  mean  term  on  the  right  is  the  average  of  the 
three  non-linear  terms  on  the  left.  We  are  in  fact  going  to  neglect  these 
non-linear  terms  but  some  justification  for  this  step  is  certainly  required. 
For  instance  the  term  ^  (uv)  is  zero  at  the  wall  eind  zero  outside  the 
sublayer,  but  reaches  a  peak  value  in  the  inner  portion  of  the  sublayer  at 
about  U^  y/v  «  10,  (See  Figure  7)*  We  can  make  this  justification  a 
posteriori  by  finding  a  solution  without  the  non-linear  terms  and  then 
comparing  the  magnitude  of  ^  (uv)  with  the  linear  terms  that  have  been 
retained.  At  this  point  we  will  state  the  results.  The  acceleration 
term  ^  is  the  leading  linear  term.  At  the  point  where  U^  y/v«10, 
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is  about  15^  of  ^  (^)  .  At  the  same  point  the  linear  convective 

terms  such  as  U(^/^)  are  about  of  It  is  evident  that  the 

non-linear  terms  are  significant  thou^  smaller  than  the  linear  terms. 
When  the  balance  of  terms  is  examined  over  the  frequency  spectrum  it  is 
found  that  at  the  low  end  of  the  frequency  spectrum,  at  say  10  c.p.s., 
th.  «lue  of  may  te  comparable  In  magnitude  to  the  terns 

that  are  retained.  But  linearization  of  the  equations  appears  to  be  a 
reasonable  first  step  towards  a  theory.  Of  course  the  non-linear  terms 
would  be  essential  in  any  theory  of  turbvilence.  But  here  the  turbulent 
field  at  the  edge  of  the  sublayer  is  assumed  to  be  known  from  the  e}q>erl- 
ments.  Our  purpose  is  merely  to  represent  the  fluctuation  field  between 
the  known  field  at  the  edge  of  the  sublayer  and  the  wall. 


niie  fluctuation  field  can  now  be  represented  by  a  superposition  of 
Fourier  components,  each  component  of  which  can  be  separately  analyzed. 
The  type  of  disturbance  that  will  be  used  will  be  based  on  our  interpre¬ 
tation  of  the  experimental  data.  As  we  have  already  suggested  one  result 
of  the  analysis  will  be  that  the  terms  containing  the  viscosity  are  only 
slgnlficeuit  in  a  narrow  region  near  the  wall.  It  is  this  region  where 
the  viscous  terms  are  important,  that  we  Identify  as  the  sublayer  of  the 
turbulent  flow.  It  shoiild  be  emphasized  that  we  are  not  concerned  with 
the  energy  balance  of  the  frequency  components  as  is  the  case  for  the 
study  of  oscillations  in  a  laminar  flow.  No  stability  calculations  are 
Involved  in  the  description  of  the  sublayer. 
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5.  SIMPLIFIED  TEIEORy 


A  field  of  turbulence  can  be  represented  by  a  Fourier  superposition 
of  elementary  plane  vorticity  waves.  The  velocity  vector  associated 
with  a  given  vector  wave  number  lies  in  a  plane  perpendicular  to  the 
wave  vector.  On  the  other  hand,  we  "have  previously  shown  that  most  of  the 
fluctuation  energy  near  the  wall  is  an  induced  field  arising  from  vorticity 
distributed  throu^out  the  boimdary  layer.  The  type  of  elementary  dis¬ 
turbance  that  is  appropriate  then  depends  on  the  nature  of  the  large  scale 
vorticity  in  the  boundary  layer,  and  not  on  the  vorticity  field  at  the 
edge  of  the  sublayer. 


Physically,  one  might  expect  the  shearing  action  of  the  mean  flow  to 
stretch  out  vortex  lines  in  the  direction  of  the  mean  flow.  There  might 
then^be  a  preference  in  the  large  scale  motions  for  vortices  neajrly  paraULel 
to  the  wall.  Since  the  root-mean- square  w'  fluctuation  is  approximately 
.7  of  the  root -mean- square  u'  fluctuation  (Figure  6),  some  type  of  "three 
dimensional"  dlsturbeuice  would  appear  to  be  necessary.  Accordingly,  we 
have  chosen  to  asstime  a  simple  form  of  oblique  disturbance  at  the  edge 
of  the  sublayer  6^  as  shown  in  Figure  8.  We  might  imagine  this  disturbance 

a 

to  be  associated  with  a  'periodic  yortlcitv  normal  to  the  d^r*f^ot-!nn 


J-T..  ^  _ »J 


_ 1 _ _ 
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rather  than  normal  to  the  I  direction  as  would  be  true  for  a  shearing  wave. 
The  velocity  q  varies  as  q  =  Q  cos  (2jf/x,e)|,  where  w/u  =  tan  0.  As  shown 
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energy  containing  eddies  is  large.  Iftider  these  conditions,  if  we  consider 
a  disturbance  of  a  particular  wave  length,  the  variation  of  the  perturbation 
velocity  with  y  in  the  vicinity  of  the  sublayer  can  be  neglected  cQmi»,red 
to  the  variation  Inside  the  sublayer. 


^Is  oblique  dlsturbemce  is  carried  downstream  with  the  velocity  U 
in  the  x  direction.  Thus  the  wave  length  in  the  x  direction  X  =  X./  cos  Q. 

X  s 

Introducing  complex  notation,  we  can  write  with  p  =  2jtf . 


C  ^(v 

I 

U  =  TRp  <n_e  > 

p  J 

I  i(k  •x-0t)'' 

w  =  Re  Bj^e  "  > 

L 


where  U^  =  >  tan  ©  =  and  Re  steind  for  Read  Part.  Slmila-rly 

the  fluctuating  pressure  field  can  he  represented  as 


f  i(k  .x-pt)  ^ 
P  =  Re  -^P^e  V 


where  p^  is  complex. 


Xuside  the  sublayer y  we  have 


with 


r  i(k  .x-pt)^ 

u  =  Re-^  h(y)e  V 

L  J 

w  =  Re  <  k(y)e  > 

L  J 

C  i(k..x-pt)1 

V  =  Re  •|g(y)e  ^  j. 

u=v=w=Oat  the  wall  y  =  0. 


Now  close  to  the  wall  where  u  -►0,  and  v -«•  0,  the  convective  terms  in  the 
equations  of  motion  ceui  be  neglected.  For  instance,  in  the  sublayer 


oU/^ox  =  ik^.u,  uu/ut  =  "  ipu,  so  thal 


11^  =  ^  • 

I  Btl 

Also  from  the  continuity  equation  with 


=  lie  u, 
X  ^ 


=  ik^.w.tajiQ,  so  that  v  =  -  j*  ik^udy  -  j  ik^wtanOdy 


or  y  =  0(  I  ®)  I  )  - 
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Provided  the  angle  0  does  not  approach  v  =  0(  I  k^.u.y  I  ). 

This  limitation  on  6  does  not  appear  to  be  significant.  We  are  interested 
in  the  u  and  v  components  at  the  wave  number  as  obtained  in  a  one¬ 
dimensional  spectral  analysis.  The  wave  nianber  for  the  oblique  disturbance 
k|  =  k^/ cos  0  so  that  k^  -•  oo  as  ©  -•  9©°  •  Since  the  energy  spectrum  for 
the  fluctixation  field  falls  off  very  rapidly  with  increasing  frequency, 
the  high  frequency  oblique  waves  where  0  -•  90°  probably  make  a  negligible 

contribution  to  the  u  and  v  fields.  With  t—  =  0  f— ),  we  have 

ay  'y" 


Since  U  ®  .8u,  for  the  large  scale  eddies  the  convective  terms  can  be 

W  X 

neglected  for  y  0.  As  we  shall  show  later,  the  thickness  of  the  viscous 
region  for  a  disturbance  of  frequency  f  is 

2  2  /  2  \ 

)  •  Then  =  0  frequency  of 

of  interest  only  the  derivatives  with  respect  to  y  need  be  retained.  The 
equations  of  motion  can  then  be  simplified  to  the  following  form 


ciy 


together  with  the  continuity  equation.  It  is 


evident  that  to  this  order  of  simplification  there  is  no  coupling  between 
the  u  and  w  fluctuation  fields  and  the  u  component  can  be  solved  for 
separately. 
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Again  under  the  restriction  that  the  disturbance  angle  9  does  not 
approach  90°,  it  can  now  be  shown  that  the  terms  in  equation  (6)  are  of 
higher  order  compared  to  the  terms  in  eqmtion  (5)»  For  instance 


=  •  i3v  =  0 


[(k^.3.|)(l  +  tan^0).y]  so  that  g  =  0  (k^.|)(l+tan^9) 

fel 


or  a  nift.Yinmtn  of  0(k  6  )  (provided  9  is  reasonably  limited)  at  the  edge  of 

X  s 

the  viscous  region. 


Now  k^6g  =  0  With  -  1.25  x  10^  cm/sec,  =  5,5  x  10“^ 


■  ^w 


at  f  =  500  c.p.s.  This  merely  expresses  the  fact  that  near  the  wall  v 
is  of  higher  order  than  u. 

An  Important  consequence  is  that  the  pressvire  field  does  not  vary 
with  y  in  the  viscous  region.  We  have, 


It  also  follows  that  ^ 

oy  ox 


•ir  (“  ■^)  =  ik  —  and  since  — 

^  'p  X  p  ^  P  oy  'ot 


±  Ss.  ,  o(^) 


(4-©  “  0(k4.p.S  (l+tan^e).,) 

oy  p  cot  A  d 


or  finally 


/I  ^ 
"S  /dUN 


o  ■ 

^  2  [l+tan^9 

=  V.  •  — 7^ - 


/) 


Similarly  it  can  be  shova  that 


1  9p 

.p 

rawH 

.  mi 

=  0  (k 


2  1 1  +  tan'^O 
c  2 


y^) 


Thus  the  pressure  field  associated  with  the  disturbance  is  constant 

2  2 

through  the  viscous  region,  since  k  5  «  1, 
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To  find  the  solution  to  equation  (5),  we  represent  the  disturbance 
as  the  sum  of  two  ccmponents,  where  u  =  u^.  The  component 

represents  the  disturbance  velocity  before  account  is  teiken  of  the  direct 
wall  effect.  As  previously  discussed,  we  assume  that  the  variation  of 
u^  with  y  in  the  vicinity  of  the  viscous  region  is  negligible  compared 
to  the  variation  of  u  that  occurs  inside  the  viscous  region.  Since  the 
convective  terms  have  been  neglected  u^  is  not  affected  by  the  mean  flow 
and  is  constant  throughout  the  viscous  region,  u^  represents  the  dis¬ 
turbance  velocity  component  directly  associated  with  the  wall  friction. 
Thus,  at  the  edge  of  the  viscous  region  u^  ~“0,  and  at  the  wall  u^  =  -  ti^. 
If  we  write  equation  (5)  as 

a(u,+  U,)  ,  ^^(u^+  u  ) 

— 5t - — 7T - 

dy 

then  the  equations  for  u^^  and  u^  throughout  the  viscous  region  are 
dt  p  dx 

^  :rT- 

cjy 


where  as  we  have  already  shown,  ^  ^  does  not  vary  through  the  viscous 
region.*  Substituting 

■i  f  \r 

j.  ^  J.V  •rk.— pu  / 

^  =  ^2  (y)e 

..  iphp 

we  2iave  h^  +  =  0  which  has  the  simple  solution 


h^  =  -  e 


SO  that  finally  u  =  Uj^  +  u^  =  Re 


■dV 

'e 


I 


This  same  separaujLon  of  equation  w&is  used  by  rrancox  in 

discussing  oscillations  in  a  laminar  boundary  layer. 


It  is  evident  that  the  region  of  rapid  chamge  in  u  is  0(^/ v/f)  as 
previously  asserted  so  that  the  extent  of  the  viscous  region  is  different 
for  each  frequency,  decreasing  in  size  as  the  frequency  is  increased. 

In  order  to  conipare  these  results  with  experiment,  we  calculate 

i?  =  ^  Re  |ji(y)  h* (y^  .  Introducing  the  dimensionless  variable 

(8)  —  =  l-2e  ^cos  Y+e 

Cl  4 


This  function  is  shown  in  Figure  9.  For  each  frequency  component 
"7  2 

u  /(Cj^  ^2^^  1  at  Y  «  5*  rapid  decrease  of  u  occurs  for  0<Y<2. 

Entering  the  viscous  region  u  /(C^  /g)  first  increases  reaching  a  peak 
value  at  Y  =  2.2.  Examination  of  the  details  of  the  solution  shows  that 
this  increase  of  disturbance  level  entering  the  vlscoxas  region  arises 
because  near  the  outer  edge  of  the  viscous  region  and  Ug  have  an  in- 
phase  component  instead  of  being  l80°  out  of  phase  as  at  the  wall. 

Figures  5  eind  9  can  now  be  used  to  calculate  the  variation  of  the 
root-mean- square  fluctuation  level  u'  near  the  wall.  According  to 
Figure  "bhe  energy  spectra  are  similar  between  y/6  =  .05  and  y/6  =  .58* 
Then  the  disturbance  energy  for  each  frequency  outside  of  the  viscous 

2 

region  is  given  by  setting  C,  /g  =  P-^  at  y/&  =  .05  from  Figure  5«  For 

u 

T 

each  frequency,  the  variation  of  u  with  y  is  obtained  from  Figure  9» 

Q3ie  u  spectra  for  various  values  of  y  can  then  be  ccanputed,  and  are  shown 
in  Figvire  10.  As  y/6  increases,  the  spectra  approach  the  spectrum  outside 
of  the  sublayer.  Integrating  for  each  value  of  y/8  and  taking  the  square 
root,  we  finally  obtain  the  variation  of  u'  near  the  wall.  In  Figure  11(a), 
the  ordinate  u'/^'^  ratio  of  u'  inside  the  viscous  region  to  the 

value  of  u'j^  outside.  According  to  the  calculations,  u'/u'^^-*  1  at  about 
y/6  =  .035»  Ibe  rise  in  the  experimental  u'  fluctuation  level  near  the 
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wall  appears  to  start  for  .05  <  y/5  <  .4,  but  it  is  not  possible  to 

establish  this  point  with  any  precision  because  of  the  variation  of  u’ 

outside  the  sublayer.  A  comparison  of  Laufer's  pipe  data  with  similar 

theoretical,  calculations  is  shown  in  Figure  11(b).  In  both  cases  the 

theory  correctly  predicts  the  total  extent  of  the  sublayer,  and  also 

the  fact  that  there  is  a  rapid  decrease  of  close  to  the  wall. 

On  the  basis  of  these  calculations,  we  feel  justified  in  identifying  the 

viscous  region  in  the  theory  with  the  sublayer  in  the  experiments.  On 

the  other  hand,  the  rise  in  fluctuation  level  entering  the  sublayer  is 

much  greater  than  the  small  rise  found  in  the  theory.  It  should  be  noted 

that  by  y/8  =  .01,  the  ratio  U./U  which  is  a  measure  of  the  relative 

w 

magnitude  of  the  convective  terms  that  have  been  neglected  is  approximately 
.6,  increasing  to  .75  by  y/b  =  .055.  Thus  the  agreement  between  theory 
and  experiment  oiig^t  to  be  better  near  the  wall. 

As  shown  in  Figure  10,  there  is  a  marked  change  in  the  si)ectral 
distribution  approaching  the  wall.  Most  of  the  energy  is  taken  from  the 
largest  scale  motions  resulting  in  a  fairly  flat  spectrum  close  to  the 
wall.  One  spectral  measurement  was  also  made  by  KLebauioff  at  y/8  =  .0011, 
deep  in  the  sublayer.  These  measurements  are  also  shown  in  Figure  10 
where  they  can  be  compared  with  the  calculated  spectrum  for  y/8  ^  .0011. 
Qualitatively,  the  measurements  confirm  the  expectations  of  the  theory. 
However,  there  is  a  significant  difference  in  magnitude  between  theory 
and  experiment  at  y/8  =  .0011.  A  somewhat  different  way  of  applying  the 
theory  suggests  itself. 

Since  we  expect  theory  to  be  better  near*  the  wall,  we  can  calculate 

the  spectra  in  the  sublayer  based  on  the  experimental  spectrum  at 

y/8  =  .0011  rather  than  the  spectnmi  outside  of  the  sublayer.  That  is, 

2 

using  Figure  9,  the  variation  of  /p  with  f  is  chosen  so  that  the 
theoretical  and  experimental  spectra  coincide  at  y/b  *  ,0011.  Then  the 
corresponding  spectra  can  be  computed  at  other  values  of  y/8.  The  results 
are  shown  in  Figure  12  and  the  corresponding  variation  of  u'/u'j^  is  shown 
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in  Figure  11(a),  Here  of  course,  u'/u'j^  for  theory  and  experiment  have 
been  set  equal  at  y/6  =  .0011.  Bie  theoretical  variation  of  u’/u'^^  is 
now  in  much  better  accord  with  the  experiments  in  the  inner  part  of  the 
sublayer,  at  least  up  to  y/6  =  .OO5.  Beyond  y/8  =  .005,  the  theoretical 
curve  rapidly  departs  from  the  experimental  data  with  the  theory  rising 
to  a  much  higher  value  of  u'  outside  the  sublayer. 


A  comparison  of  the  experimentaJL  spectrum  at  y/6  =  .05  and  the 
theoretical,  spectrum  at  y/6  =  .OO5  (Figure  12)  indicates  that  the  higher 
fluctijatlon  level  at  y/6  =  .OO5  is  due  to  an  increase  in  the  energy  in 
the  large  scale  eddies.  This  is  consistent  with  the  theory  in  that  at 
y/6  =  .005,  all  frequencies  above  f  =  6OO  c.p.s.  should  still  be  outside 
of  the  viscous  region.  Therefore  it  seems  probable  that  the  spectra 
shown  in  Figure  12  are  a  better  representation  of  the  spectra  close  to  the 
wall  than  the  spectra  in  Figure  10. 

T 

While  it  is  possible  to  determine  the  variation  of  u  in  the  viscous 
region  without  specifying  the  way  the  disturbance  energy  varies  with  the 
angle  0,  a  knowledge  of  the  three  dimensional  character  of  the  field  is 
necessary  to  determine  v  .  Suppose  we  consider  the  disturbance  velocities 
Uq  ^d  v_  associated  with  an  oblique  disturbance  at  angle  9.  Since 
equations  (5)  and  (7)  are  identicsU.  in  form  and  “  ^  =  (p  tanO,  the 

relation  w_  =  u^  tanO  holds  throughout  the  viscous  region.  We  can  then 
readily  find 
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=  [ 


1  + 


2 

tan  9 


i]  Re  I 


g(y) 


i(k^.x-St) 


} 


SO  that 

(9) 

where  T(y) 


T(Y) 

-  sinY)-2Y+2Y^ 
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Here  is  the  magnitude  of  at  the  edge  of  the  viscous  region. 

Since  the  problem  of  properly  accounting  for  the  distribution  of 

the  oblique  disturbanees  will  not  be  considered  in  this  paper,  We  are 

2 

not  able  to  compute  v~  from  the  theory. 

However,  some  remarks  about  the  application  of  the  present  theory 

to  the  determination  of  the  field  can  be  made.  The  factor  f  in  the 
2  2 

equation  for  v  shifts  the  v  spectrum  to  higher  frequencies  as  compared 

with  the  u  spectrvm.  This  filtering  action  is  a  fundamental  aspect  of 

~2 

the  response  of  the  sublayer  to  the  u  fluctuation  field.  Whereas  the 
one  dimensional  u^  fluctmtion  field  decreases  monotonically  from  f  =  0, 
the  V  spectrum  at  the  edge  of  the  sublayer  shoiild  exhibit  a  maximum, 
liiat  is,  the  v  field  close  to  the  wall  may  have  a  preferred  wave  length. 

There  are  two  further  difficulties  in  the  theoretical  determination 
of  V  which  should  be  mentioned.  In  order  to  canpute  the  total  v  ,  the 
Integration  woiild  have  to  be  carried  out  to  frequencies  greater  than 
500  c.p.s.  This  csuinot  be  done  in  a  satisfactory  way  for  the  following 
reasons.  First  the  variation  of  with  frequency  at  the  higher  frequencies 
is  not  known.  Second  a  basic  limitation  is  imposed  by  the  form  of  dis- 
txarbance  assumed  in  the  theory.  We  have  assumed  that  for  a  given  frequency 

+  V»  ^  •TP>r»  r'i.  i  T  ^  » ■*•1  «  w-A  ^  ^  ^ 
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viscous  region  can  be  neglected  compared  to  the  variation  inside  the  very 
small  viscous  region.  Since  the  thickness  of  the  viscous  region  6^  is 
ffiven  bv  6  =  0^  /I  \  thp  n  nf*  t.hp  1r  -hVipn 

r  'u  -'v  f '  -  - - s - 

.  At  f  =  500  c.p.s.,  ^  « 150.  Above  500  c.p.s.,  U..  must  necessarily 


decrease  below  .8u^,  but  there  are  insiifficient  measurements  to  establish 
the  variation  of  U„  with  f.  (This  point  will  be  mentioned  again  in 
discussing  the  wall  pressure  f luctiaations . )  Suppose  we  adopt  the  artificial 
picture  of  an  eddy  of  scade  L,  travelling  along  the  wall  with  the  average 
velocity  of  its  center  located  at  y  =  l/2.  Then  for  f  =  5.5x10^, 

L/6g  «  25,  This  guess  suggests  that  while  L/6^  will  decrease  with 
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increasing  frequency,  it  may  still  be  large  enou^-»  However,  at  say 
y/6  =  .OOS,  for  a  frequency  of  5.5x10'  ,  the  distance  frcan  the  waJ.1  y 
is  of  the  same  order  as  the  eddy  scale  L  and  it  is  obvious  that  the 
variation  of  the  perturbation  velocity  outside  of  the  viscous  region 
woiold  have  to  be  taken  into  accoimt.  Accordingly,  the  theory  may  have  to 
be  limited  in  its  application  to  the  v  spectra  below  f  =  500  c.p.s. 

We  aJ-so  find  that  the  determination  of  the  shear  stress  depends  on 
a  knowledge  of  the  distribution  of  oblique  disturbances.  According  to 
the  simplified  theoiy 
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We  recall  that 
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then  the  shear  stress  at  the  edge 


dv„ 

and  rr-  +  tt—  =  0,  where  we  have  set  v,  =  v^  =  0,  at  y  =  0.  At  the 
ooc  oy  -  “ 

edge  of  the  viscous  region  for  any  frequency,  u^  —0,  and  v^  -•  constant. 
If  we  designate  the  constant  by  V2  ^ , 

of  the  viscous  region  is  due  to  the  coupling  of  u^  and  ^200 
function  S(Y)  is  shown  in  Figure  15  together  with 

/  I  o  \ 

V  u^  /  (c,72)  . 
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It  is  evident  that  in  the  inner  portion  of  the  viscous  region  the 
relative  variation  of  u'  and  uv  with  y  for  each  frequency  ccmponent  is 
similar  to  the  overall  variation  of  u ’  and  uv  shown  in  Figure  7 • 

Increase  in  the  shear  stress  lags  behind  the  Increase  in  the  fluctuation 
level.  In  both  experiment  and  theory,  at  the  peak  of  the  u'  fluctuation 


uv  is  .7  *8  of  the  shear  stress  outside  the  viscous  region.  However, 

the  magnitude  of  uv  given  by  the  simplified  theory  is  definitely  too 
low.  This  can  be  shown  in  the  following  way. 

Sufficient  data  axe  given  by  KLebanoff  to  calculate  the  spectral 
variation  of  uv/u'v'  at  y/6  =  .05.  The  results  of  the  calculations  are 
shown  in  Figure  14.  For  comparison,  the  variation  of  uv/u’v'  inside 
the  viscous  region  according  to  the  simplified  theory  is  shown  in  Figure 
15 •  Apparently  the  correlation  coefficient  does  not  reach  high  enough 
values  in  the  viscous  region  at  least  at  the  lower  frequencies.  !Ilie 
convective  terms  that  have  been  neglected  in  the  simplified  theory  ought 
to  play  a  significant  role  in  controlling  the  phase  angle  between  u  and 
V.  Therefore  this  deficiency  of  the  theory  is  not  surprising. 

Again  the  factor  shifts  the  uv  spectrimi  to  hi^er  frequencies  as 
compared  with  the  u  spectrum  but  not  as  much  as  for  the  v  spectrum. 
Consequently  the  hl^  frequency  limitations  on  the  theory  that  have  already 
been  discussed  should  not  be  so  serious.  But  in  any  case,  a  computation 
of  uv  would  have  to  account  for  the  three  dimensional  nature  of  the  dis- 
ttirbance  field. 
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6.  MICROSCALES 


Important  changes  in  the  dissipation  derivatives  or  microscales 
occur  on  entering  the  sublayer.  Here  we  will  consider  two  of  these 
derivatives  which  can  be  discussed,  at  least  to  a  limited  extent,  using 


The  longitudinal  microscale  can  be  written  as 
(^.k/|Re{h(y).h*{y)}. 


the  simplified  theory 

(H) 


While  there  are  no  data  in  general  for  the  variation  of  the  dls- 

turbemce  velocity  with  f  above  f  =  500  c.p.s.  a  first  guess  for  this 

variation  can  be  made  at  y/5  =  .05.  For  separation  distances  S  .5  cm 

the  autocorrelation  and  longitudinal  correlation  curves  of  KLebanoff 

appear  to  coincide.  This  would  suggest  that  for  L  =  .5  cm.  or 

f  «  2x10^  c.p.s.  Previoiasly,  we  have  set  an  upper  limit  of  f  =  500 

c.p.s.  for  =  .8U^.  We  assme  then  that  the  disttirbance  velocity  varies 

linearly  from  .8u,  at  f  =  500  to  U  =  U  at  2000  c.p.s.  Above  f  =*  2000 

c.p.s.,  U  is  set  eaual  to  U..  Integrating,  we  find  (^)*'  =  5*8 

W  *  £  TT  -  CJX 

^1 

which  ccxapares  with  a  value  of  5»5  from  the  experiment.  Deep  in  the 
sublayer  the  experimental  value  for  (■^y  increases  by  about  a  factor 
of  2,  but  the  uncertainty  about  the  variation  of  with  f  is  too  great 
to  proceed  with  any  theoretical  calciilations  in  this  region. 


Thp»  microscale  transverse  to  the  wall  can  be  written  as 


=  i  Re 

'dy'  ^ 


'fh'(y).h'  (y'J’ 


(y).h  (y)[  which  becomes 
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I  — I  -  i— jr-j  e  w  ^  Outside  the  sublayer,  at 

^dy''  V  '  2 

y/S  a  .05,  we  have  assumed  a  form  of  disturbance  where  ^  =  0,  where  in 

reality  ^  is  Bmal l  but  not^zero  outside  the  viscous  region.  The 

-hi  2  ^ 

experimental  variation  of  - ^  shown  in  Figure  I6.  Approaching 


the  wall,  there  is  a  very  rapid  increase  in  the  gradient  normal  to  the 
wall.  A  rapid  increase  in  (■^)  is  also  predicted  by  the  theory.  In 
fact,  according  to  equation  (12),  (^:)  will  he  a  maxlmiim  at  the  vail. 
Calculations  have  been  carried  out  at  several  points  in  the  inner  part 
of  the  sublayer  based  on  the  measiired  spectrum  at  y/8  =  .0011,  smd  the 


jk  w  wxjLcsbv  uxxc^  ojau^a.j.x  J.cru. 
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unaerestitnate  s  the  value  of  (^)  at  y/5  =  .005,  the  inner  limit  of  the 
experimental  measurements.  A  possible  reason  for  this  discrepancy  can  be 
suggested.  The  large  gradients  in  u  In  the  simplified  theory  are  confined 
to  the  region  very  close  to  the  vail.  Eb^erimentally,  as  shown  in  Figure 
11,  substantial  gradients  in  u  are  found  over  a  much  greater  extent  of 
the  sublayer.  We  might  anticipate  that  the  simplified  theory  wovild 
underestimate  (^)  away  from  the  wall. 

The  much  greater  experimental  value  of  (^)  as  compared  with 
fir.e.e  not  rsflsct  s.  distortion  of  the  small  scale  eddies. 


/_  ^  _  -  An  *  P 

at  y/6  =  .005,  the  theoretical  spectrum  for  shows  that  practically 

all  the  contribution  to  (-s:;;)  comes  from  frequencies  <  300  c.p.s. 
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7.  THE  PRESSURE  FIELD 


Extensive  measurements  of  the  pressure  fluctuations  at  the  boundaries 
of  turbulent  flows  have  been  made  by  Willmarth  (1958>1959)*  He  found  that 
the  major  contribution  to  the  pressure  fluctixations  comes  from  large  scale 
fluctuations.  Of  particular  interest,  his  space-time  correlation  measvire- 
ments  show  that  the  pressure  pattern  is  convected  downstrejun  with  a  speed 
of  .82U^.  This  observation  is  in  good  eigreement  with  the  present  theory 
where  the  large  scale  flucttiations  move  downstream  at  the  mean  velocity 
of  the  middle  region  of  the  boundary  layer.  In  fact  measurements  of  the 
wall  pressure  field  may  provide  a  means  of  experimentally  establishing 
the  variation  of  disturbance  velocity  with  frequency  for  the  hi^er  fre¬ 
quencies  in  the  sublayer.  Corcos  and  Winkle  (i960)  have  found  that  by 
mafr-trig  a  spectreuL  resolution  of  the  longitudinal  space-time  correlation 
a  functional  relationship  between  convective  velocity  and  frequency  can 
be  fo\md.  As  we  would  expect  theoretically  the  higher  frequencies  are 
convected  more  slowly  than  the  low  frequencies.  However,  theoretical  cal- 
ciilatlon  of  the  spectrum  emd  the  magnitude  of  the  pressure  fluctiiations  at 
the  boundary  layer  is  another  matter. 


The  pressure  fluctuations  in  a  turbulent  shear  flow  may  be  much  larger 
than  in  a  field  of  Isotropic  turbiolence  at  comparable  tvirbulent  fluctuation 
levels.  This  is  the  case  near  the  edge  of  the  laminar  sublayer  in  a  turbu¬ 
lent  boundary  layer.  Very  large  pressure  fluctuations  are  associated  with 
the  linear  terms  in  the  equations  of  motion.  Just  outside  of  the  sublayer, 
retaining  all  the  linear  terms,  equation  (l)  becomes 
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The  fluctuation  field  at  the  edge  of  the  sublayer  is  now  represented  by  a 
superposition  of  Fourier  components.  Using  the  same  notation  as  in  the 
previous  sections,  we  can  write  equation  (15)  as 


i(k^x-3t) 
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Then  the  spectral  variation  of  p  /p  is  given  by 


(14)  ^  =  F 
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In  isotropic  turbulence,  U  =  U.,  =  0,  and  the  pressure  field  due  to 

~  dy 

these  linear  terms  vanishes.  At  high  frequencies,  U  =  U,  aad  the  first 
and  third  terms  do  not  contribute  to  the  pressure  field.  However,  most 
of  the  contribution  to  the  pressure  fluctuations  ccmies  from  the  low  fre¬ 
quencies  and  these  terms  are  Important.  The  spectral  fimctions  F  auad 
P  -2  have  been  measured  at  the  edge  of  the  sublayer  at  y/8  =  .05,  ^ 

v  2 

and  can  be  directly  used  in  calcvilatlng  p  /p  .  An  additional  asstmntlon 


inus’t  be  made  In  orde'T*  detennine  F  -j—  . 

”  (iu)v 

For  a  given  frequency  component,  the  shear  correlation  coefficient 

uv/u'V  =  cos  0  where  0  is  the  phase  angle  between  the  u  and  v  velocities. 

If  we  represent  u  on  the  positive  read,  eocis  in  the  complex  plane,  then  a 

negative  correlation  coefficient  indicates  that  v  is  either  in  the  second 

or  third  quadrants  since  0  must  lie  between  90°  and  270°.  Now  according 

to  the  simplified  theory,  l80°  ^  0  <.270°  throughout  the  viscous  region 

(measiiring  0  coimterclockwise  from  u).  Assuming  then  that  v  lies  in  the 

third  quadrant  the  correlatiaa  spectrum  F  , .  »  can  be  computed  using  the 

( lu^v 

variation  of  the  experimental  phase  angle  0  with  frequency  as  given  in 
Figure  14.  That  is  if  y  is  the  phase  angle  between  iu  and  v,  then 
COS7  =  sin  0,  and  is  negative  making  the  third  term  in  eqviation  (l4) 

positive.  The  resxilting  spectrum  for  p  /p  is  shown  in  Figure  17  where  we 
have  set  =  .Su^^.  If  v  lies  in  the  second  Instead  of  the  third  quadrant, 
then  the  sign  of  P  is  reversed  although  the  magnitude  is  the  same. 

In  that  case  the  pressure  fluctviation  spectrum  would  be  represented  by 
the  lower  boundary  of  the  cross  hatched  region  in  Figure  17.  The  second 
teim  of  equation  14,  v  /k  (dU-/dy)"  is  also  given  in  Figure  17  to 
indicate  what  might  be  expected  from  the  linear  theory  at  the  higher 
frequencies.  • 


It  is  interesting  to  compare  these  calculations  vith  the  pressure 
spectnan  for  an  isotropic  field,  where  the  pressure  fluctuaticais  are  due 
to  the  non-linear  terms.  Batchelor  (1956,  p.  ifla.)  has  given  an  ex¬ 
pression  for  the  pressure  fluctuations  in  an  isotropic  field  in  terms 
of  the  three  dimensional  energy  spectrum.  The  first  step  then  is  to 
find  an  appropriate  energy  spectrum  for  this  flow.  Since  the  Reynolds 
number  of  the  turbulence  is  large,  the  one  dimensional  spectrum  F 
at  v/B  =  .0*5  can  be  fitted  by  the  function  ^ 
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where  L  is  the  longitudinal  integral  scale  (Dryden  19^3)*  Then  using 

X 

the  transformation  between  the  one  and  three  dimensional  energy  spectra 
for  isotropic  turbulence,  the  three  dimensional  spectrum  P(k)  can  be 
computed  using  Batchelor’s  theory  where 
00 


I 


2 

P(k)  dk  =  p  /p  .  Finally,  the  one  dimensional  pressure 


2 

spectrtm  P<w,  where  we  represent  p  /p  by  jT,  can  be  obtained  from 


Cl^  -  I 

T  J 


00 


p(k) 


dk  and  is  also  shown  in  Figure  17 


(transformed  to  a  frequency  spectrum).  The  values  given  are  appropriate 
to  an  isotropic  field  with  each  component  equal  to  u'.  It  is  evident 
that  at  the  edge  of  the  sublayer,  the  linear  terms  are  the  main  source  of 
the  pressiire  field. 


The  present  theory  provides  a  clear  basis  for  having  a  pressxire  field 
at  the  boundary.  For  each  frequency  component  ^(^)  «  0  across  the 

1  ckirc&iM  'CViron  +-V^/^nrr‘h  -f*!  nr»+.no  *♦- "f  nn  orirri‘nri'ntfari‘hfl  vnn'f  aVi  n.*h  flVlia  WP.11* 


the  pressure  field  does  not  vanish  with  them.  This  is  in  contrast  to 
an  isotropic  pressure  field  where  p  /p  OC  (u*)  •  It  would  then  follow 
that  the  press\jre  fluctuations  at  the  boundary  should  be  of  the  same 
order  as  the  pressure  fluctiiations  at  the  edge  of  the  sublayer.  Using 
boundary  layer  parameters  appropriate  to  Klebanoff's  experiment,  Willmarth’ 
wall  luBasurBiuBxi'ts  arB  also  shown  in  FiguTB  IT •  InuBgraoing  "tiiB  spBctra 


in  Figure  17,  and  taking  the  square  root,  we  have 


■n  f 

r _ 


pu. 


=  50,  for  the 


linear  terms,  3.7  for  the  isotropic  field  and  2.4  from  Willmarth.  It  is 


apparent  that  the  wall  pressure  field  is  much  smaller  than  what  would  be 
e3q)ected  frcm  the  present  theory.  Perhaps  the  non-linear  terms  that  have 
been  neglected  are  essential  for  the  computation  of  the  pressure  fluctu¬ 
ations  in  the  sublayer.  In  part  this  may  be  because  a  major  contribution 
to  the  pressure  field  comes  frcm  the  lowest  frequencies  and  at  these 
frequencies  the  non-linear  terms  are  relatively  more  important. 


8 ,  TURBULENCE  PRODUCTION 


A  good  desJ.  of  significance  has  been  attached  to  the  experimental 
observation  that  the  flow  of  turbulent  energy  from  the  mean  flow  to  the 
fluctuating  field  reaches  a  maximum  in  the  sublayer.  This  observation 
has  sioggested  the  following  physical  picture.  Adjacent  to  the  wall,  there 
is  a  laminar  flow.  Moving  out  from  the  wall  a  region  of  instability 
develops  because  of  the  high  disturbance  level  of  the  turbul  ence  in  the 
boundary  layer.  Along  the  outer  edge  of  the  "laminar  flow”,  transition 
to  tvirbiilence  occiirs,  producing  a  region  of  intense  turbulence  and  of 
course  '’tui’bulence  production".  Finally  going  still  further  from  the  wall, 
the  flow  settles  down  to  a  reasonably  well  behaved  fully  tiarbulent  flow. 


In  terms  of  the  theory  given  in  this  paper,  the  peak  in  turbulence 
production  does  not  suggest  such  a  physical  model.  Furthermore  if  we 
start  with  the  condition  that  the  total  shear  stress  for  the  mean  flow 


be  constant  across  the  sublayer,  then  we  can  readily  show  that  there  will 
be  a  peak  in  turbulence  production  where 


—  dU 
-  p  uv  =  p  ^ 


The  total  shear  is  -  p  uv 


dU  ^  4.  4. 

^  ==  =  T  =  Constant, 
dy  w 


The  turbulence  production  (Pr)  is  given  by 

X  —  dU  >  /dUx  dU 

(Pr)=-p  uv  •5—  =  (t  -  p  (— )  — 
^  ^  dv  '  V  ^  '  d  V  dv 

%/  --»/  - V 


Then 


^  o  dU  /d^Ux 

(Pr)  =  -  2p  ^  (-^)  = 

dy  *'  dy 


0, 


and  we  have  p  ^  for  maximum  production.  This  result  is  confirmed 

by  the  experiments  of  KLebanoff  and  Laufer.  Thus  while  the  region  of 
maximum  tvirbulence  production  is  important  for  determining  the  structure 
of  the  turbulence,  it  need  not  be  given  any  phenomenological  significance. 


Generally  speaking  in  the  fully  developed  part  of  the  boundary  layer, 
the  energy  flows  frcsa  the  mean  flow  into  the  large  scale  eddies.  At 
y/6  =  .05,  from  KLebanoff,  70^  of  the  shear  stress  is  found  between 
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0  <  f  <  500,  which  compares  with  80^  of  u  below  f  =  500.  This 

dU  w 

situation  is  altered  in  the  sublayer.  At  the  point  where  1^-  ^ 
the  loss  in  turbulent  shear  stress  has  come  from  the  large  scale  eddies. 
This  means  that  the  turbulence  production  goes  into  ever  smaller  eddies 
as  the  wall  is  approached.  Therefore  the  region  of  maximum  turbulence 
production  may  not  be  important  in  determining  the  energy  balance  of 
the  large  scale,  energy  containing  eddies  in  the  boundary  layer. 
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9.  lAMINAK-TUREULENT  TRANSITION  IN  STOONG  TINfflUIMCE 

The  most  satisfactory  correlation  of  the  exi)erimentaLl  data  for 
the  Reynolds  number  of  boundary  layer  transition  in  strong  free  stream 
turbulence  has  been  given  by  Taylor  (1956).  Taylor  derived  his  transition 
parameter  on  the  assumption  that  the  fluctuating  pressiore  gradients  of  the 
turbulence  cause  momentary  separation  of  the  laminar  boimdary  layer,  thereby 
leading  to  tremsitlon.  In  recent  years,  detailed  investigations  of 
transition  have  been  carried  out  by  Schubauer  and  KLebanoff  (1955).  In 
no  case  do  they  find  any  experimental  support  for  the  ides  that  momentary 
separation  is  involved  in  the  transition  phenomena. 

The  flow  near  the  wall  of  a  laminar  boundary  layer  with  a  strong 
free  stream  turbvilence  Is  in  msiny  respects  similar  to  the  flow  in  the 
sublayer  of  a  turbulent  flow.  The  scale  of  the  free  stream  turbulence  is 
in  general  large  compared  with  the  thickness  of  the  laminar  boundary  layer. 
The  free  stream  turbulence  moves  downstream  with  the  free  stream  flow,  and 
therefore  at  a  disturbance  velocity  much  larger  than  the  local  mean 
velocities  in  the  boundary  layer  near  the  wall.  In  accordance  with  the 
present  theory,  we  should  eijqpect  to  find  a  "sublayer"  of  the  eddies  of 
the  free  stream  turbulence  in  a  smsill  region  close  to  the  wall.  Ilhe 
neglected  non-linear  terms,  which  are  of  concern  in  the  treatment  of  the 
sublayer  of  a  turbulent  flow,  are  of  no  importance  in  this  case. 

The  present  theory  has  been  applied  to  the  experimental  conditions 
given  by  Dryden  (1956).  These  data  are  at  a  transition  Reynolds  number 
based  on  the  displacement  thickness  6*  of  R^^  w  500,  which  appears  to  be 
the  lowest  transition  Reynolds  number  to  be  foiand  in  the  literature,  ttiis 
is  very  close  to  the  minimum  critical  Reynolds  number  of  Rj.„  =  400  of  the 
small  disturbance  stability  curve  for  laminar  boundary  layer  oscillations. 
Accordingly,  it  should  not  be  necessary  to  take  into  account  the  amplifi¬ 
cation  of  the  disturbance  level  in  the  boundary  layer  in  this  case. 
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Llepmarm,  Lavifer  and  Liepmann  (1951)  have  measured  the  energy- 
spectrum  for  the  free  stream  turbiilence  in  a  vind  tunnel  under  experi¬ 
mental  conditions  reasonably  close  to  those  of  Dryden,  We  estimate 
L^/u^  «  1x10  ^  for  Dryden  where  =  .9x10  ^  for  Liepmann.  Hius 

Liepmann's  spectrum  can  be  used  to  calculate  the  "sublayer"  for  Dryden 's 
experimen-bal  point.  It  is  found  that  the  disturbances  due  to  the  free 
stream  turbulence  are  damped  down  by  -viscosity  in  a  "sublayer"  approximately 
.15  cm  thick.  The  close  similarity  be-bween  the  physical  picture  for  the 
sublayer  of  a  turbulent  flow  and  the  "sublayer"  of  the  free  stream  tur¬ 
bulence  in  a  laminar  boundary  layer  is  shown  in  Flgijre  l8.  The  scale 

—2 

corresponding  to  50^  of  the  u  energy  is  shown  in  each  diagram.  It  seems 
clear  that  the  physical  model  used  by  Taylor,  in  which  the  mean  boundary 
layer  responds  to  the  pressure  gradients  associated  wi-fch  the  turbulence 
is  not  applicable  to  the  description  of  the  laminar  boundary  layer  in  a 
strong  free  stream  tiirbulence. 

Nevertheless,  Taylor's  parameter  does  appear  to  correlate  -the  a-vailable 
experimental  da-ba.  We  woiild  now  like  to  show  that  by  introducing  an 
additional  physical  assumption  we  can  ob-bain  a  transition  parameter  very 
similar  to  that  of  Taylor's.  No  attempt  will  be  made  to  justify  this 
assumption  here  al-bhough  it  has  led  to  reasonable  results  in  some  pre-vlous 
work  (Lin,  1955,  P.90) 

We  assume  that  -bhe  onset  of  transition  depends  on  the  relative 
amplitude  of  the  Reynolds  stress  associated  wi-th  the  txirbulent  fluctuations 
£ind  the  shear  in  the  mean  flow.  This  is,  we  adopt  as  a  rough  criterion 
for  the  onset  of  transition,  -the  requirement  that 
_  , ^  f  \  .  .  W 

-  M  uv  =  j.  •iJ-  g 

Prom  the  simplified  theory,  the  shear  stress  outside  the  "sublayer",  is 
of  the  form 

-puvCt  p(u')^ 

w 
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Here  U  -  U  >  the  free  stream  velocity.  Ihe  frequency  f  can  be  replaced 
w  1 

by  f  =  where  represents  the  "scale"  of  the  turbulence.  Then 


-  p  uv  OC  p(u'  ) 


.^2  ^ 


Now  the  thickness  of  the  laminar  boundary  layer  at  transition  is 


/x,  V 

UX-  ^  U 


Substituting,  we  have 


P( 

Finally,  we  obtain 

.u' 


u'f  .  ^  (Re,p.  ^ 


’i  / 

7. 


V 


Vvx 


tr 


X.  '  4 


x  e 


Taylor's  paj:ameter  is 


V. 


(-1^)  =  f-^Re.  )  where  L  is  the  lateral  Integral 

'‘U  '  'L  '  tr'  y 

1  y 


scale  of  the  free  stream  turbulence.  The  available  transition  data  are  too 
scattered  to  make  it  possible  to  distinguish  between  a  l/4  or  l/5  power 
variation  with  x^^L.  From  the  point  of  view  of  correlating  the  experimental 
data  either  parameter  would  be  equally  effective  although  Taylor's  parameter 
is  based  on  the  assumption  that  transition  is  caused  by  the  small  eddies, 
whereas  for  the  new  parameter  it  is  assumed  that  transition  is  caused  by 
the  larse  eddies. 


10.  CONCLUDING  REIMARKS 


The  vorticity  field  responsible  for  the  tiirbulent  fluctuations  in 
a  boundary  layer  or  similar  shear  flow  is  swept  along  with  the  velocity 
of  the  fluid  elements.  The  velocity  fluctuation  field  associated  with 
this  vorticity  field  is  altered  by  the  wall  in  two  ways.  GThe  effect  of 
the  boTjndary  condition  v  =  0  is  to  increase  the  magnitude  of  the  wall 
velocity  fluctuations  u  and  w  in  the  plane  of  the  wall.  The  induced 
velocity  at  the  wall  associated  with  each  element  of  vorticity  is  doubled 
by  the  image  vortex  element  required  to  cancel  v  at  the  wall.  This  form 
of  "wall  effect"  extends  across  the  boundary  layer  and  beyond  into  the 

potential  flow.  As  a  result  of  the  boundary  conditions  u  =  w  =  0  the 

turbulent  velocity  fluctuations  are  directly  damped  down  by  viscosity  in 
a  thin  layer,  the  sublayer. 

The  equations  of  motion  for  the  tiirbulent  velocity  and  pressure 
fluctuations  are  applied  only  in  this  narrow  viscous  region.  The  aim  of 

the  theory  is  to  say,  in  detail,  how  a  known  turbulent  field  dies  out  to 

nothing  at  the  wall. 

A  simplified  form  of  the  theory  is  given  in  this  paper.  Only  the 
leading  terms  in  the  differential  equations  are  retained.  Furthermore 
while  the  three  dimensional  character  of  the  fluctuation  field  is  recog¬ 
nized  and  introduced  at  any  early  stage  the  calculations  are  not  carried 
far  eno\jgh  to  make  the  three  dimensionality  important. 

The  same  basic  approach  can  be  used  to  develop  a  more  accurate 
description  of  the  sublayer.  The  linear  convective  terms  can  be  retained 
in  the  equations  and  solutions  are  being  obtained  using  computing  machines. 
Also  the  Influence  of  the  distribution  of  oblique  disturbances  on  the 
fluctuation  field  is  being  explored.  It  is  hoped  to  report  the  results 
of  these  studies  in  the  future. 

Calculations  with  the  theory  could  be  applied  to  such  qviantities  as  v 
and  uv  if  the  variation  of  disturbance  velocity  with  frequency  at  the 
higher  frequencies  could  be  established.  This  might  be  done  from  space- 
time  correlation  measurements  of  wall  pressure  fluctuations. 

kl 


Finally  there  are  related  problems  in  vhlch  the  present  aK>roach 
may  be  useful.  We  have  already  shown  that  the  theory  applies  to  the 
"sublayer"  of  the  free  stream  turbulence  in  a  laminar  bomdary  layer. 

A  similar  approach  might  be  useful  in  the  description  of  turbulent  heat 


and  mass  transfer. 
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y/8  =  .58  From  Klebonoff 
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FigTJre  1.  Typical  Boundsiry  Layer  Ener^^y  Spectrum 
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Figure  2.  Effect  of  Frequency  Cutoff  and  Disturbance  Velocity  on  Avito- 
correlation  Coefficient. 
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Figtire  3,  Comparison  of  Autocorrejlation  suad  Longitudinal  Correlation 
Coefficients  From  Fa-'/re. 
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Figure  6.  Variatioin  of  Fluctuation  Levels  Across  Boundary  Layer. 
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Fluctuation  Level  and  Turbulent  Shear  Stress  Near  the  Wsill 
'ftaom  Klebanoff  and  Laufer. 
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Figure  7* 


tz 


Figitre  8.  Assumed  Form  of  Oblique  Distxirbance  at  Edge  of  Sublayer. 
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Figure  10. 


Energy  Spectra  In  the  Sublayer  Based 


on  Spec'fcmjn  Outside  Sublayer* 
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Figure  11.  (a)  Compsirison  of  The' 

in  the  Sublayer  - 
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Theory  using  spectrum  at  y/8=.0011 
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and  ExperiLmejnt  for  Turbulent  Fluctuations 
.ndary  Layer. 
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Figure  u.  (b)  Comparison  of  Theory  and  Experiment  for  Ttirbulent  Fluctioations 
in  the  Sublayer  -  Pipe  Flow. 
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Figure  12.  Energy  Spectra  in  the  Sublayer  Based  on  Spectrum  Inside  Sublayer. 
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Figijre;  l4. 


Variation  of  Shear  Correlation  Coefficient  with  Frequency  Outside 
sifter  KLebanoff. 
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Figxore  15 .  Theoretical  Variation  of  Shear  Correlation  Coefficient  in  the 
ViscoiiB  Region. 
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Figure  Comparison  of  G?heory  and  Experiment  for  Triansverse  Microsctile  in 

the  Sub'Layer. 
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The;  so-ca]LLed  lamlntar  sublayer  is  shiown  to  be  the  region  where  the  turbuLenl:, 

'  velocity  fluctuations  are  dJ.rectly  dissipated  by  vi.scosity,  A  sjjnplified 
.  I:lnearl2.ed  forca  of  the  equations  of  moticn  for  ithe  tvirbulent  fluctuiitlons  is  used 
to  describe  th<»  turbulent  f  j.eld  be'twe(!n  the  wall  and  the  fully  turbulent  part  of 
tlie  flow\  'Eie  mean  flow  In  the  sublayer  and  th<;  turbulence  field  outside  the  sub¬ 
layer  are  assumed  to  be  jtno^m  from  the  e3?perimeiats .  The  thlcikness  of  the  subilayer 
arrises  naturally  in  the  theory.  It  is  shown  that  tihe  large  sicale  fiLuctuatlons 
containing  iinosi:  of  the  turbulent  enerfjy  a, re  convected  downsti-eam  with  a  velocity 
:  diaracteristic  of  the  middle*  of  the  bcaundary  layer.  Thus  Tailor’s  ]::ypotheslE 
'  does  not.  apply  to  these  ;iarge  scale  fluctuations  near  the  waJl.  CaJLculatlons  are 
g:Lven  fcr  the  <margy  spectrei  and  u’  f]Luctuation  le'vel  in  the  sutlay(?r  and  other 
aispects  of  the  fluctuation  field  are  cLlscussed.  Eaanlnlng  ttie  effect  of  strong 
I  f:ree  str-eam  tuiibulence  on  iBuninar  boundary  layeir  ti-ansitlon,  it  appiars  that  the 
;  pliyslcal  model  underlyings  Teiylor’s  pajiameter  is  incorrect. 


Tlie  e;o- called  Jamlnar  sul)layer  is  shcAsn  to  be  the  region  ^^here  the  tuit>ulent  j 
velocity  fluctuations  are  dlreictly  dissipated  by  ^kriscoslty.  A  simplified  i 

linearized,  form  of  the  equations  of  motion  for  the  turbulent  f]LuctuationiB  is  used  | 
to  descrihe  the  turl)ulent  field  between  the  wall  (and  the  fully  turbulent  paa-t  of 
the  flow.  The  m<ian  flow  in  the  sublayer  and  the  turbulence  field  outsldie  the  sub-  ! 
layer  fire  assumed  to  be  known  from  the  experlmentis .  The  thiclaieBS  of  the  8\iblayer  ; 
arises  nat.urally  in  the  theorj^.  It  is  shovrn  that  the  large  scale  fluctuations  1 
containing;  most  of  the  turbuleint  energy  are;  convected  doivnst retun  with  a  'irelocity 
characteristic  of  the  middle  of  the  boimdaiy  layer.  Thua  Taylor's  hypothesis  | 

does  n(?t  a.pply  tiD  ttiese  large  scale  fluctus.tions  laear  the  \rall,  CalciilatioiiB  are 
given  jfor  the  enijrgjr  spectra  tind  u'  fluctus.tion  level  in.  the  siiblayer  and  other  ! 

aspects  of  the  fluctuation  field  Eire  discusaed.  lExaminlng  the  effect  of  strong  1 

free  stream  turbiolerice  on  laialLnai*  boundary  layer  transition,  lit  appears  that  the 
physictil  model  ujiderlylng  Taylor's  peirameter  is  Ijacorrect.  | 
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Tlie  Eio-call<ed  3.amlnar  sublayer  is  shovfn  to  be  the  region  where  the  turliulent 
velocity  fluctuations  are  directly  dissipated  by  ^viscosity.  A  simplified 
linearized,  form  of  Wie  equations  of  motion  for  the  turbuilent  fluctuations  i(5  used 
to  descrile  the  tur^)ulent  field  between  the;  wall  land  the  fully  turbulent  part  of  i 
the  flow.  Tbe  m(“an  flow  in  tlie  sublayer  and  the  turbulence  fi<5ld  outside  tlie  sub-  I 
layer  are  assumed  to  be  known  from  the  experiments.  The  thickness  of  the  sublayer  i 
arises  nat;urally  in  the  theor^f^.  It  is  shovn  that  the  large  scale  flu.ctuatlons  1 
containing;  most  of  the  turbul<;nt  energy  are?  convected  dcwnstreim  with,  a  velocity 
characteristic  of  the  middle  of  the  boundaiy  layer.  ThiJis  Taylor’s  hypothesis 
does  not  apply  to  tliese  large  sca.le  fluctusitlons  near  ttie  wall.  Calculations  are 
given  for  the  energ^r  sp-ectra  and  u'  fluctuaitlon  level  in  the  sublayer’  and  either 
aspects  of  the  fluctuation  fl<5ld  are  discufised.  :Exaa1nlng  the  effect  of  strong 
free  stream  turb^ilence  on  lam^Lnai*  boundary  layer  transition,  it  appears  that  the 
physlail  model  uadej'lylng  Taylor's  parametasr  is  Incorrect. 


